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Traditionally, phase transitions are defined in the thermodynamic limit only. We 
discuss how phase transitions of first order (with phase separation and surface ten- 
sion), continuous transitions and (multi)-critical points can be seen and classified 
for small systems. "Small" systems are systems where the linear dimension is of the 
characteristic range of the interaction between the particles. I.e. also astrophysical 
systems are "small" in this sense. Boltzmann defines the entropy as the logarithm of 
the area W{E,N) — e'^'^'^' of the surface in the mechanical N-body phase space 
at total energy E. The topology of S{E, N) or more precisely, of the curvature 
determinant D{E, N) = d'^S/dE'^ * d^S/dN^ - {d^S/dEdNf allows the classifica- 
tion of phase transitions without taking the thermodynamic limit. Micro-canonical 
thermo-statistics and phase transitions will be discussed here for a system coupled 
by short range forces in another situation where entropy is not extensive. The first 
calculation of the entire entropy surface S{E, N) for the diluted Potts model (ordi- 
nary {q = 3)-Potts model plus vacancies) on a 50 * 50 square lattice is shown. The 
regions in {E, N} where D > correspond to pure phases, ordered resp. disordered, 
and D < represent transitions of first order with phase separation and "surface 
tension". These regions are bordered by a line with D = 0. A line of continuous 
transitions starts at the critical point of the ordinary {q = 3)-Potts model and runs 
down to a branching point Pm- Along this line VD vanishes in the direction of 
the eigenvector vi of D with the largest eigen-value Ai « 0. It characterizes a 
maximum of the largest eigenvalue Ai. This corresponds to a critical line where the 
transition is continuous and the surface tension disappears. Here the neighboring 
phases are indistinguishable. The region where two or more lines with D = cross 
is the region of the (multi)-critical point. The micro-canonical ensemble allows to 
put these phenomena entirely on the level of mechanics. 

PACS numbers: 05.20.Gg, 05.50-fq, 05.70Fh 



I. INTRODUCTION 

Nuclei, atomic clusters and astrophysical objects are not large compared to the range of their 
forces. Therefore, these systems are not extensive. (In the following we call systems "non- 
extensive" in a somewhat more general sense: If they are divided into pieces, their entropy is 
not the sum of the entropies of their parts in contrast to conventional extensive systems where this 
is assumed at least if the pieces are themselves macroscopic). Although the largest systems possible 
belong to this group we call these systems "small" . But also at phase transitions in systems with 
short range forces does the entropy of the surface separating the different phases not scale with 
the size of the system. Lateron we define systems to be "small" precisely by the condition that 
the entropy S{E, N, V) does not scale with the number of particles or the volume and especially if 
S{E, N, V) has some convex regions. 

Most applications of thermodynamics to "small" systems are more or less transcriptions from the 
thermodynamics of macroscopic systems, c.f. the book by Hill Conventional thermo-statistics, 
however, relies heavily on the use of the thermodynamic limit (V^oo\pf/v, or ^ const.) a-nd extensiv- 
ity, c.f. e.g. the book of Pathria |^]. This is certainly not allowed for our systems. Extensivity is 
nowadays considered to be an essential condition for thermodynamics to work, c.f. Lieb and Ygn- 
vason That the micro-canonical statistics works well also for "small" systems without invoking 
extensivity will be demonstrated here for finite normal systems. The use of the thermodynamic 
limit and of extensivity, however, is closely intervowen with the development of thermodynam- 
ics and statistical mechanics since its beginning more than hundred years ago. When we extend 
thermodynamics to "small" systems we should establish the formalism of thermodynamics start- 
ing from mechanics in order to remain on a firm basis. This is an old program since Boltzmann 



and Einstein [g|j5|. We will see how this idea guides us to more and deeper insight into the most 
dramatic phenomena of thermodynamics, phase transitions. Moreover, it gives the most natural 
extension of thermo-statistics to some non-extensive systems without invoking any modification of 
the entropy like that proposed by Tsallis |^ see also . 

In the following section we sketch a deduction of thermo-statistics from the principles of mechan- 
ics alone. Nothing outside of mechanics must be invoked. This was the starting point of Boltzmann 
, Gibbs Hi, Einstein and the Ehrenfests [|lO|,0 at the beginning of this century. They 
all agreed on the logical hierarchy of the micro-canonical as the most fundamental ensemble from 
which the canonical, and grand-canonical ensembles can be deduced under certain conditions. Ac- 
cording to Gibbs the latter two approximate the micro ensemble in the thermodynamic limit of 
infinitely many particles if the system is homogeneous. Then surface effects and fluctuations can be 
ignored relatively to the bulk mean values. This is the main reason why the thermodynamic limit 
became basic in the statistical foundation of macroscopic thermodynamics. However, it was Gibbs 
who stressed that the equivalence of the three ensembles is not even true at phase transitions 
of first order, even in the thermodynamic limit. 

This chapter repeats the familiar deduction of statistical mechanics as it was first formulated 
by Gibbs ||. However, we will indicate at which point the main stream of logical steps towards a 
thermodynamics of "large" systems has to be left in order to cover also "small" and non-extensive 
systems. 

Just a remark is neccessary here: One might think that a small system embedded in a heat bath 
may be described by the canonical or even the grand-canonical ensemble. A heat bath interacts 
with the system via its surface. For small systems these surface mechanisms are important. They 
affect the physical properties and must be specified. Certainly more must be controlled in such 
cases than just the mean energy (temperature) as is assumed in a canonical treatment. This would 
be only correct for an infinite system, with short range interactions, and if we wait an infinitely 
long time, then the details of the physical mechanisms in the surface contacting the system with 
the bath become unimportant relatively to the bulk. 

The third section addresses phase transitions. First, the basics of the conventional definition of 
phase transitions in the thermodynamic limit (^^— >oo|n/v, or /x const.) by the theory of Yang and 
Lee are reminded. Then this definition is translated into the micro-canonical picture and the 
source of the non-analyticities of the canonical thermodynamic potentials is traced to the topology 
of the micro-canonical entropy function, esp. its curvature. 

In the fourth section we explain the general features of the micro-canonical phase diagram as 
function of fixed energy density (e — E/V) and particle number- in — A^/l/)/magnetization- 
density, at first, with the well known example of the Ising model and the [q = 1) diluted Potts 
model. It is shown how its phase space is limited by the line of random configurations which have 
the maximum possible entropy. Further the region of "field driven" phase separation into sponta- 
neously magnetized clusters (first order transition) can be well studied in the energy-magnetization 
or the energy-particle number phase diagram. (Remark: the distinction between "field driven" and 
"temperature driven" transition makes sense only for the Ising model. Due to its intrinsic symme- 
try M — > —M the energy- axis is also the main curvature direction for M = 0. This is of course 
not the case for systems without this symmetry.) Then the diluted Potts model with q = Z is 
introduced. It is sketched which features of its phase diagram are to be expected in the case of a 
"small" system. 

Section V gives a short discussion of the main steps of the simulation and in the section VI the 
numerical results are shown. We present the first view of the entire entropy-density surface s(e, n) 
as function of energy density (e) and particle number density (n) for the diluted (g = 3)-Potts 
model on a finite 2-dim lattice with "volume" y = More insight into the phase diagram of 
the model is given by an overall view of the determinant of the curvature density d{e,n). Various 
kinds of phase transitions and critical lines can be seen and classified. A subsection is addressed 



^In the following we often skip the volume V as third variable as this is always hold fixed. Of course for 
finite systems the shape of the container matters. In the diluted Potts model we use a square lattice with 
periodic boundaries ("square torus"). 
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to the essential differences between the micro-canonical and the grand-canonical phase diagram. 
The information lost by the latter representation is emphasized. 

The following section discusses the validity of the second law of thermodynamics in non-extensive 
systems with a convex s(e,n). The consequences of convexity of the entropy for Wcinhold's ge- 
ometric interpretation of thermodynamics are discussed. It will be explained that a convexity of 
S{E, N) is not in contradiction to the second law of thermodynamics. 

Finally, our main conclusions on phase transitions in "small" non-extensive systems are summa- 
rized in section VIII. 

II. ON THE MECHANICAL BACKGROUND OF THERMO-STATISTICS 

An equilibrated many-body system is characterized by few macroscopic quantities: 

1. Its energy E, mass (number of atoms) N, volume V, 

2. its entropy S, 

3. its temperature T, pressure P, and chemical potential fi. 

There are important qualitative differences between these three groups: All variables of the first 
group have a clear mechanical significance. They arc conserved and well defined at each point 
of the N-body phase space. The internal dynamics of the system cannot leave the shell in phase 
space which is defined by these variables. Also entropy as the most important quantity within 
thermodynamics has a clear mechanical foundation since Boltzmann. His gravestone has the 
famous epitaph: 



relating the entropy S to the size W{E, N, V) = eotrS{E — Hn) of the energy (E) surface in the 

N-body phase space at given volume (V). Here eo is a suitable small energy constant which does 
not affect any variation of the entropy, is the A^-particle Hamiltonian, and 



The set of points on this surface defines the micro-canonical ensemble. In contrast to the conserved 
quantities which are defined at each phase space point, the entropy refers to the whole micro- 
canonical ensemble. 

Remark: For a system with discrete energies Ei e.g. a lattice or a quantum system with energy 
spacing eo one should define the micro-canonical partition sum by the number of states at this 
energy. When we discuss derivatives of W we imagine a suitable smoothing of this. 

It is important to notice that Boltzmann's and also Einstein's formulation allows for defining the 
entropy entirely within mechanics by Smicro '■= ln[W{E, N, V)]. It is a single valued, non-singular, 
in the classical case miiltiply differentiablc, function of all "extensive" , conserved dynamical vari- 
ables. No thermodynamic limit must be invoked and this definition applies to non-extensive like 
our "small" systems as well. 

The third group of quantities which characterize the thcrmodynamical state of an equilibrated 
many-body system, temperature, pressure and chemical potential have no immediate mechanical 
significance. Within micro- canonical statistics they are defined by the derivatives of the entropy 
S{E, N, V) vs. the conserved quantities : 



S=k- InW 




(1) 



}__dS_ 
T~dE 



(2) 




(3) 



(4) 
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From the mechanical point of view they are secondary, derived quantities. This difference to the 
two other groups of variables will turn out to be significant for "small" systems. Again, like entropy 
itself, these quantities characterize the whole micro-canonical ensemble, not an individual point in 
the N-body phase space. 

Starting from this point, the conventional thermo-statistics assumes extensivity and explores 
the thermodynamic limit (1^—>oo|n/v, or ^ const.) c.f. [Q. This procedure follows Gibbs ||^. He 
introduced the canonical ensemble, which since then is the basic of all modern thermo-statistics. 
The link between both ensembles is established by a Laplace transform. E.g. the usual grand- 
canonical partition sum is the double Laplace transform of the micro-canonical partition sum 

Z(r, V) = ir — dN e-yE-^N~TSiE,N,v)yT ^ Yl fr de dn e-^[-pn-Ts(e,„,y)]/T_ 
i/o ^0 Co i/o 

In the thermodynamic limit it is useful to work with the energy density e = E/V , the particle 
number density n = N/V , and the entropy density s — S/V, cf. the right expression in cq.(|5|). 



III. PHASE TRANSITIONS MICRO-CANONICALLY 

According to Yang and Lee phase transitions are indicated by singularities of the grand- 
canonical potentials (cx yln[Z]) as function of z = e^/"^ on the positive real z-axis. These, however, 
can occur in the thermodynamic limit (V^^ 00 1^ const.) only. For finite volumes the number of 
particles N is finite. Consequently, Z is a sum of a finite number of powers and ^ln[Z] is 
analytical for positive z at any T. Are there no phase transitions in finite systems? 

There are phenomena observed in finite systems which are typical for phase transitions. Some- 
times this is even so in astonishingly small systems like nuclei and atomic clusters of ~ 100 atoms 
-|l6|. In ref. we showed that their characteristic parameters as transition temperature, la- 
tent heat, and surface tension are - in the case of some metals - already for thousand atoms close, 
though of course not equal, to their known bulk values. Therefore, it seems to be fully justified to 
speak in these cases of phase transitions of first order. 

We need an extension of thermodynamics to "small" systems which avoids the thermodynamic 
limit. However, here is a severe problem. The three popular ensembles, the micro-canonical, the 
canonical, and the grand- canonical ensembles are not equivalent for "small" systems. The energy 
per particle can fluctuate around its mean value <E/N> in the (grand-)canonical ensemble whereas 
the energy fluctuations are zero in the micro-canonical ensemble. Moreover, the heat capacity is 
strictly positive in the canonical ensembles whereas it may become negative in the micro ensemble. 
It was Gibbs himself who warned for the use of the canonical ensemble at phase transitions of first 
order In this situation it is certainly advisable to keep close contact with mechanics. It is 

helpful to realize that the fundamental micro-canonical ensemble as introduced by Boltzmann is 
the only one which has a clear mechanical definition ||l8|,|| for finite systems. 

To extend the definition of phase transitions by Yang and Lee to finite systems we study 
which feature of the micro-canonical partition sum W{E, N, V) leads to singularities of the grand- 
canonical potentials yln[Z] as function of z = e^/"^ by the Laplace transform cq.(||). In the 
thermodynamic limit F — > cxd|^ const, this integral can be evaluated by asymptotic methods. As 
far as the entropy surface s(e,rt) has everywhere negative curvatures the integrand of eq. (H) has 
a single maximum. For large V the Laplace integral (r.h.s of eq.(^)) is then dominated by the con- 
tribution of this peak. This is a stationary point {es,ns} where = ds/de, v — —[ijT = ds/dn. 
If there is only a single stationary point then there is a one to one mapping of the grand-canonical 
ensemble to the micro-canonical one and energy-fluctuations disappear. 

This, however, is not the case at phase transitions of first order. Here the grand-canonical 
ensemble contains several Gibbs states (stationary points) at the same temperature and chemical 
potential which contribute similarly to the integral (||). Consequently, the statistical fluctuations 
of e and n do not disappear in the grand-canonical ensemble even in the thermodynamic limit. 
This is the reason why Gibbs himself excluded phase separations Between the stationary 

points s(e,n) has at least one principal curvature > 0. Here van Hove's concavity condition |19| 
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for the entropy s{e,n) is violated. In the thermodynamic hmit these points get jumped over by 
the integral and ln[Z] becomes non-analytic. Consequently, we define phase transitions also 
for finite systems topologically by the points and regions of non-negative curvature of the entropy 
surface s(e,n) as a function of the mechanical, conserved "extensive" quantities like energy, mass, 
angular momentum etc.. 

The central quantity of our further discussion the determinant of the curvatures of s{e,n) is 
defined as 



d{e, n) 



de^ dnde 
dedn dn-^ 



(6) 



Also critical fluctuations, i.e. abnormally large fluctuations of some extensive variable in the 
grand-canonical ensemble or the eventual divergence of some susceptibilities are micro-canonically 
connected to the vanishing of the curvature determinant, e.g. in the following examples of d{e,n) 
or d(e, m) respectively: 



d 



The micro-canonical specific heat is given by : 

de _ Snn 

dT ^ ^ T^d(e,n) ' 

d{f3iy) 

d{en) 

or the isothermal magnetic susceptibility by : 
dm 

as 



with s 



J. d{e,m) ' 



ee 



dede 



etc. 



(7) 
(8) 

(9) 
(10) 



In the case of a classical continuous system s(e, n) is everywhere finite and multiply differentiable. 
In that case the inverse susceptibilities [cmicro{e,n,V)]~^ and [Xmicro,T{G,n,V)]~^ are well be- 
haved smooth functions of their arguments even at phase transitions. Problems arise only if the 
susceptibilities are considered as functions of the "intensive" variables T, and v or B [EO . In 
the case of lattice systems we can only assume that the inve rse susceptibilities are similarly well 
behaved. This will be further illuminated in the section VI B. 



Experimentally one identifies phase transitions of first order of course not by the non-analyticities 
of yln[Z] but by the interfaces separating coexisting phases, e.g. liquid and gas, i.e. by the 
inhomogeneities of the system which become suppressed in the thermodynamic limit in the grand- 
canonical ensemble. This fact was early realized by Gibbs and he emphasized that using 5* vs. 
volume at phase separation "has a substantial advantage over any other method because it shows 
the region of simultaneous coexistence of the vapor, liquid, and solid phases of a substance, a region 
which reduces to a point in the more usual pressure-temperature plane." That is also the reason 
why for the grand-canonical ensemble the more mathematical definition of phase transitions JlSf is 
needed. The main advantage of the micro-canonical ensemble is that it allows for inhomogeneities 
as well and we can keep much closer to the experimental criteria. 

Interfaces have three opposing effects on the entropy: 

• An entropic gain by putting a part (A^i) of the system from the majority phase (e.g. solid) 
into the minority phase (bubbles, e.g. gas) with a higher entropy per particle. However, 
this has to be paid by additional energy AE to break the bonds in the "gas" -phase. As 
both effects are proportional to the number of particles A^i being converted, this part of the 
entropy rises linearly with the additional energy. 

• With rising size of the bubbles their surfaces grow. This is connected to an entropic loss 
(surface entropy) proportional to the interface area due to additional correlations between 

2/3 

the particles at the interface. As the number of surface atoms is cx N-^ this is not linear in 



AE and leads to a convex intruder in S{E, N, V), the origin of surface tension |22 
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• An additional mixing entropy for distributing the iVi-particles in various ways over the bub- 
bles. 

At a (multi-) critical point two (or more) phases become indistinguishable because the interface 
entropy (surface tension) disappears. 



IV. THE MODEL AND ITS QUALITATIVE FEATURES 

In this paper we want to discuss how critical and also multi-critical points manifest themselves 
in the micro-canonical statistics of a finite diluted Potts model. 

It is helpful first to remind the main physical effects to be expected: We start with the well known 
properties of the Ising model (written here as a Potts model with q — 2 ||2^). The Hamiltonian 
and the magnetization are defined as: 

n.n. pairs 



H=- J2 Sa.:'r,-Bj2<^^ (11) 

i,j i 

M = Y,a,. (12) 

i 

The double sum runs over nearest neighbor spin-pairs only [ p3| , p4[ . 

The two-dimensional Ising model was extensively discussed as function of temperature T and 
magnetic field i? = in the grand-canonical ensemble for L ^ oo (thermodynamic limit) p^,p[. 
As known from Onsager's solution it has a continuous (second order) critical point at _B = 
and (3c ~ jr ^ 0.51n(\/2 + 1) j2^j2^. Above the critical point, T > Tc there is a disordered 
phase with vanishing magnetization for B = 0. Below Tc there is a bifurcation of the spontaneous 
magnetization curve M{T, B ^ 0) — sign{B)L^mspo7itan{T). As function of the extensive variables 
{E,M} between M — zLL'^nispontaniT) there is a whole region with phase separation which is 
inaccessible to the grand-canonical ensemble. Also the internal energy per lattice point e{T, B = 0) 
is analytically known in the thermodynamic limit c.f. This may be inverted to get the function: 

-(e,i?-0) = ( . ,e>e.^-(l + ^) 

I sign{B)mspontan[e) ,e<ec 

Then the phase diagram in the {e, m} plane can be drawn as in fig.|l|: 
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FIG. 1. Phase diagram in magnetization (m) and energy (e) per spin particle for the Ising model 
(ordinary {q — 2)-Potts model) in the thermodynamic limit. The upper parabolic line is the locus of the 
completely random configuration e — — (m^ + 1) with maximum entropy. The physical possible phase 
space points are in the two shaded regions. The darkly shaded region is the phase space of first order 
transition and coexistence of drops with positive and drops with negative magnetization. It is inaccessible 
in the grand-canonical ensemble. The dot gives the critical point (cc — ~(1 + ■^), i^c ~ 0). 

At energies e>ec = — (1 + -^) one has a disordered ("gas") phase. Atec = — (l + -^),r7i==0 
there is a critical point of continuous transition with vanishing surface tension. From here an 
inhomogeneous region of condensed phase starts towards lower energies (several clusters of positive 
or negative magnetization coexist separated by interfaces) . In the thermodynamic limit the overall 
magnetic field < B > vanishes here. This region is bordered by the curve nispontanie, B — 0). 
Approaching this border from inside with fixed e and rising m, the clusters with the opposite 
magnetization get depleted and the surfaces between the clusters disappears. (However, even 
when there are no clusters of the wrong magnetization there should still be fluctuations with 
the other, wrong magnetization because these states have a nonvanishing entropy. There is a 
general impossibility to distinguish a phase separation from a fluctuation within an individual 
configuration. This is a realistic example for the definition of a phase as a property of the entire 
ensemble only.) This region of phase separation is darkly shaded in fig.|l| and is inaccessible to the 
grand-canonical ensemble. 

The following Gedanken-experiment may illuminate why this region cannot be accessed in the 
grand-canonical ensemble: Suppose one prepares the system at a sharp energy and magnetization 
in the region with positive curvature. Then one puts the system into a heat bath with temperature 
Ttr and magnetization (particle-number in the case of the diluted Potts model) bath with ^tr 
where it will be suddenly exposed to energy- and magnetization- (particle number) fluctuations 
of arbitrary size. Because of the positive curvature of its entropy the system is unstable and 
will receive or loose sufficient energy to convert entirely into the gas or entirely into the solid 
phase where its free energg is lower. The latent heat which is neccessary for this will be spended 
(absorbed) by the bath. 

The Ising model can be modified by the transformation Oi = -^(cj + 1) and Hig = 
~5 Si j ^oi,-i-i^oj.+i into the diluted {q = l)-Potts model. Here the lattice points with 0^ = 
are vacancies. Figure (^) shows the phase diagram of the diluted {q = l)-Potts model correspond- 
ing to the phase diagram of the Ising model as discussed above. The three lines indicate (from 
left to right) the ground state {eo{n) = — 2n), the critical line, and the line of random (maximum 
entropy) configuration {emax{n) — — 2n^). The dot gives the critical point which corresponds to 
the critical end point of the phase separation in the liquid-gas transition. Between the ground 
state and the critical line is the {e,n}-region (darkly shaded) of first order phase transition with 
phase separation. Here the system prefers inhomogeneous configurations with droplets enbedded 
in the gas. This corresponds to the above mentioned similar inhomogeneous region of the Ising 
model at energies below the Ising critical point (ec = — (1 + "^)j™c = 0), where clusters of positive 
magnetization and negative magnetization coexist. In the thermodynamic limit we expect the 
micro-canonical entropy s(e, n) to be flat in the shaded region. The whole darkly shaded region of 
inhomogeneities is not accessible in the grand-canonical ensemble. This is a striking example of the 
dramatic loss of information when the grand-canonical ensemble is addressed to phase transitions. 
Above the critical line we have the region of the pure disordered phase. Here s(e, n) is concave. 
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FIG. 2. Phase diagram for the diluted {q — l)-Potts model in the thermodynamic limit. With the 
transformation at = 2oi — 1 of the two possible occupation numbers to the spin on site i and m = y^^. aj 
this is analogous to the energy-magnetization plot of the Ising model. The dot gives the critical point 
(fie « —0.853, Uc = 0.5). The darkly shaded region is the phase space of phase coexistence which is 
inaccessible in the grand-canonical ensemble. 

In the following we investigate the 3-states diluted Potts model now on a finite 2-dim (here 
= 50^) lattice with periodic boundaries in order to minimize effects of the external surfaces of 
the system. The model is defined by the Hamiltonian: 

n.n. pairs 
H = - ^ OiOjScr^^aj 
i,j 

i 

Each lattice site i is either occupied by a particle with spin ai — 1,2, or 3 or empty (vacancy). 
The sum is over neighboring lattice sites i, j, and the occupation numbers are: 

J 1 , spin particle in site i 
* \ J vacancy in site i 

This model is an extension of the ordinary {q — 3)-Potts model to allow also for vacancies. 
At zero concentration of vacancies {n = 1), the system has a continuous phase transition at 
Be = 1 + ~ 1-58 @,||. With rising number of vacancies the probability to find a pair of 
particles at neighboring sites with the same spin orientation decreases. The inclusion of vacancies 
has the effect of an increasing effective qeff > 3. This results in an increase of the critical energy 
of the continuous phase transition with decreasing n and provides a line of continuous transition, 
which is supposed to terminate when qeff becomes larger than 4. Here the transition becomes first 
order. 

At smaller energies the system is in one of three ordered phases (spins predominantly parallel in 
one of the three possible directions). We call this the "solid" phase. This scenario gets full support 
by our numerical findings in the next sections. 

In the following, we will show for the case of the diluted {q = 3)-Potts model how the total 
micro-canonical entropy surface S{E, N) uncovers first order phase transitions, continuous phase 
transitions, critical, and multi critical points even for small systems and non-extensive systems. 
This is an extension of previous work on systems with one thermodynamic degree of freedom like 
the fragmentation phase transition in highly excited nuclei , in atomic clusters ||l7|] , and in the 
ordinary Potts model |22y27f| . The main purpose for the study of systems with two thermodynamic 
degrees of freedom like {E, N} is the possibility to localize (multi-) critical points. 
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(14) 



(15) 



V. SIMULATION METHOD 



The simulation methods proceeded in two steps: We first covered aU space {E — e* L^, N — n* 
L^}, L = 50 by a mesh with about 1000 knots with distances of Ae — 0.08 and An — 0.04. At each 
knot {ci, Uk} we performed micro-canonical simulations (« 2* 10^ events) to get a histogram of the 
probabilities P(ei, rife) = e'^ for the system to be distributed in the narrow region ± 4) * ( A^^ ± 4) 
of phase space. Local derivatives (3 = {dS{E,N)/dE)j^, — /3/i = v = {dS{E,N)/dN)^ in each 
histogram give a first "measurement" of the "intensive" quantities. These are used to interpolate 
f3, V over the gaps between the knots of our mesh c.f. pst . 

1. From the first interpolation we get an estimate /3o, for the center of any region {Ae * An}. 

2. In a second iteration step we sample Pi{e,n) = P{e,n) * e~^°^ (e-Mo") ^ broader {e,n} 
domain. This is chosen to touch or overlap the next neighboring region. This way a dense 
covering of the entire {e,n} space is achieved. 

3. A fourth order regression fit fi{e, n) in e and in n to ln[Pi(e, n)] is performed, from which we 
determine the two principal axes of curvature Ai, A2 of s(e, n) and then can rotate if desired 
the frame of the basic rectangle in {e, n} to be parallel to the main curvature axes. 

4. The last step is then a longer and extensive sampling of P2{e,n) — e^^'^^'^'"''* Pi{e,n) in the 
new basic rectangle. 

5. Eventually one has to recycle steps 3 and 4 a few times up to the time one obtains a sufficiently 
uniform covering of ^2(6, n). 

6. At the end on gets back the original by P(e, n) = e-^''('='")e'^«^'(''~''o")F2(e, n) 

We obtain a broad and uniform covering of the whole {e, n} space. This way one gets the most 
effective (and fast) "measurement" of s(e,n), /3(e, n), v^e^n) up to all third order derivatives of 
s(e,n) e.g. See(ej?^) or See„(e,n) etc.. Technical details of our method will be published in [ p9[ . 
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VI. RESULTS 




FIG. 3. Entropy Smicro{e,n). 



Figure (||) shows Smicroi^, n). Grid lines are in the direction [e — eQ(n)\/[emax{n)~eo{n)\ — const 
and n =const, (eo(ri) — —2n, e„jQ^(n) = ^^^)- The grey levels are determined in fig.^: Regions 

above CPmB : concave, d > 0, pure phase (disordered, gas), in the triangle APmC concave, pure 
phase, ordered (solid); Below AP,nB: convex, d < 0, phase-separation, first order; At the dark 
lines like APmB we have d(e,n) ~ 0: these are termination lines of the first order transition; 
Medium dark lines like CPm'- vi • Vd = 0, here the curvature determinant has an extremum in the 
direction of the largest curvature eigenvector vi (remember: in the normal concave region of s(e, n) 
the curvatures are both negative. A vanishing or even small positive curvature is an abnormal large 
curvature) ; In the cross- region (light gray without grid) we have: d = A Vd = this is the locus 
of the multi-critical point Pm where s(e,n) is (numerically) flat up to at least third order in Ae 
and An. 

The convex region of phase separation corresponds to the similar region in the Ising lattice gas, 
respectively the original Ising model as function of magnetization, c.f. the darkly shaded region 
in figure (||). The main new feature of the phase diagram of the diluted {q = 3)-Potts model 
compared to the diluted {q = l)-Potts model is the critical line from C to Pm- 

In figure (^) the determinant of curvatures of s{e,n): 
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is shown. On the diagonal we have the ground-state of the 2-dim Potts lattice-gas with cq — — 2n, 
the upper- right end is the complete random configuration (here without contour lines), with the 
maximum allowed excitation e^and = region above the line CPmB we have the 

disordered, "gas". Here the entropy s{e,n) is concave {d > 0), both curvatures are negative (we 
have always the smaller one A2 < 0). This is also the case inside the triangle AP,nC (ordered, 
"solid" phase). In these regions the Laplace integral eq.(^) has a single stationary point. They 
correspond to pure phases. 



11 



FIG. 4. Contour plot of the determinant of curvatures d{e,n) defined in eq.(|l^).The grey levels are as 
in fig.^ Regions above CPmB : concave, d > 0, pure phase (disordered, gas), in the triangle APmC 
concave, pure phase (ordered, solid); below APmB: convex, d < 0, phase-separation, first order; At the 
dark lines APmB we have d(e, n) = 0: termination lines of the first order transition; Medium dark lines 
e.g CPm.: vi • Vd = 0; here the curvature determinant has a minimum in the direction of the largest 
curvature eigenvector Vi; in the cross-region (light gray) we have: d = 0AVd = O this is the locus of the 
multi-critical point Pm where s(e,n) is fiat up to at least third order in Ae and An. The two horizontal 
lines give the positions of the two cuts shown in figsJaJq. 
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FIG. 5. Cut through the determinant d(e, n) 
along the line shown in figure (^)^at const. 
n = 0.95, through the critical line CPm close to 
the critical point C of the ordinary Potts model 
(n ~ 1) 



FIG. 6. Cut through the determinant d(e, n) 
along the line shown in figure (^ at const. 
n = 0.57, slightly below the multi-critical region. 
There are several zero points of the determinant of 
curvatures: The left one is simultaneously a max- 
imum with Vd = and consequently critical as 
discussed above 

Below APjnB s{e,n) is convex (d < 0) c.f. figure (0), corresponding to phase-separation, first 
order. At these {e,n} the Laplace integral (^) has no stationary point. Here we have a separation 
into coexisting phases, e.g. solid and gas. Due to the intra-phase surface tension or the negative 
contribution to the entropy by the additional correlations at the phase boundaries (surface), s(e, n) 
has a convex intruder with positive largest curvature. In ||2^;|l^ it is shown that the depth of the 
convex intruder in s(e,n) gives the surface tension. 

At the dark lines like APmB we have d{e,n) = 0. These are the termination lines of the first 
order transition. At these lines one of the two phases is depleted and beyond all particles are in 
the other phase (solid or gas respectively). 

Along the medium dark lines like PmC we have vi • Vc? = 0, here the curvature determinant 
has a minimum in the direction of the largest curvature eigenvector Vi. The line PmC towards 
the critical point of the ordinary {q = 3)-Potts model at e = —1.58, n — 1 correponds to a critical 
line of second order transition which terminates at the multicritical "point" Pm- It is a deep valley 
in d{e, n) c.f. fig|| which rises slightly up towards C. On the level of the present simulation we 
cannot decide whether this rise is due to our still finite, though otherwise sufficient, precision or 
is a general feature of finite size. (The largest curvature Ai of s(e,n) has a local maximum with 

Ai ^ 0, or d ~ 0). Because of our finite interpolation width of Ae ~ ±0.04, An ^ ±0.02 it might 
be that this valley of d(e, n) gets a little bit filled up from its sides and the minimum is rounded, c.f. 
figj^. The valley converts below the crossing point P„j into a flat ridge inside the convex intruder 
of the first order lattice-gas transition see also fig.^. 

In the cross-region (light gray in figj^) we have: d = A Vd = 0. This is the locus of the multi- 
critical point Pm where s(e,n) is (numerically) flat up to at least third order in both directions Ae 
and An. It is at Cm ~ —1, rim ^ 0.6 or = 1.48 ± 0.03, Vm = 2.67 ± 0.02. Naturally, Pm spans 
a much broader region in {e,n} than in j/}, remember s(e,n) is flat near Pm- This situation 
reminds very much the well known phase diagram of a '^He -^He mixture in temperature vs. mole 
fraction of '^He c.f. fig. 3. in ref. [0. 
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Figure (0) gives a 3D- view of the entropy-surface s(e, n) — f3o{e — Cq) — lyoin — no) along a broad 
strip with constant (3 = —1.85 (through the region of phase separation, roughly parallel to the 
ground state eo(n) = —2n). The overall convex intruder of the first order transition is well seen. 
In its middle the narrow ridge where s(e, n) is again concave can be imagined. This is the origin 
of the medium dark line PmD in figures (H) and (^. 




FIG. 7. Entropy surface s' — s(e,n) — 1.85(e — eo) — 3.61(7i — no) along a strip with constant /3 = 1.85 
(roughly parallel to the ground state eo(n) — —2n through the region of phase separation). 



A. On the topology of curvatures 

The two eigenvalues of the curvature matrix (16) are: 



Al,2 - ± ^ ViSee + Snn? " 4d (17) 

and the corresponding eigenvectors are : 



1 



en 



(19) 



At critical points the following conditions hold: 



d = y^^^L-D = , (20) 
oyen) 
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Ai 


= 


A2 






=0 





Here the directions /3 =const. and /i = const, are parallel 
and we have : 





1 



va<o 



1 



(21) 
(22) 

(23) 

(24) 

(25) 
(26) 

(27) 
(28) 



The vanishing of d is not sufficient for criticahty. Physically, it means that the surface entropy 
(tension) and with it the interface separating coexistent phases disappears. This, however, can 
also signalize a depletion of one of the two phases. At a critical end-point, however, the interface 
disappears at a non vanishing number of atoms in each of the two phases. I.e. in an infinitesimal 
neighborhood of a critical point, d must remain zero. In a topologically formulation a critical 
end-point of first order transition is at: 



d = 
and 
Vi-Vd = 0. 



(29) 
(30) 



This is a generalization of the well known condition for a continuous transition in one dimension: 
the simultaneous vanishing of f3'{e) — and of the curvature of /3(e), /3"(e) — 0. 




-2 -1.5 -1 -0.5 



e 

FIG. 8. Direction of the largest principal curvature vi . 

Figure ^ shows a map of some trajectories which follow the eigen-vector vi with the largest 
curvature eigen- value Ai. In the region of the convex intruder (Ai > 0) i.e. the region of phase- 
separation Vi is ~ parallel to the ground state e — ~2n. Also the lines of P =const. and /x =const. 
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follow approximately this direction. Their Jacobian d{j3v) / d{en) = d{e,n) is negative but small. 
This reminds of the situation in the thermodynamic limit where this region of phase coexistence 
is flat, both intensive variables are constant and the Jacobian d ^ — 0. One can also see in fig.^ 
how the direction of the largest curvature Vi turns into the e-direction when one approaches the 
critical point C of the ordinary {q ~ 3)-Potts model at n = 1. 

Inspection of fig.§ shows that along the line CPm of second order transition vi • Vd = 0, i.e. 
d{e, n) has a deep and sharp valley. This line is the locus of a minimum of d in the direction 
of the largest eigenvalue Vi of the curvature. In the direction of CPm the determinant c?(e, n) is 
slightly growing towards the ordinary Potts critical point C . Figure shows an original narrow 
histogram of s(e, n), (Ae = 0.08 and An — 0.04) on the line CPm near to the point C which shows 
that the curvature of s(e, n) is rather small in the direction of the larger curvature eigenvector Vi. 

At n = 1 we know that for an infinite system the ordinary in — \) three state Potts model 
has a second order transition at e — —1.58 where the curvature of s(e) vanishes, See — 0. I.e. 
the component vi • Vc? of Vc? indicates nicely the locus of the second order "temperature driven" 
transition of the ordinary Potts model. 



N-NO 




E-EO 



FIG. 9. Histogram of s(e,n) — (e — eo)/3(eo,no) — (n — no)v{eo,nQ) at eo = —1.48, no = 0.94 in the 
rectangle eo ± 0.04, no ± 0.02 viewed in the direction of the valley PmC of d{e, n), the second-order critical 
line towards the critical point C of the ordinary q = 3-Potts model at e = —1.58. The approximate 
vanishing of the largest curvature in the direction vi (pointing from left to right) at this point eo,no is 
clearly seen. 
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B. The information lost in the grand-canonical ensemble 



1.25 




nu 3 . 5 

FIG. 10. Plot of the entropy Smicro{l3, v) as function of the "intensive" variables {v — —Pfi) in the figure 
labeled as "nu" and /3 is called "beta". The lines which build the surface are lines for /3 =const. The two 
bold ones indicate the cuts shown in figsjn] and |l|. The positions of the points A,D,B,C defined in figure 
(^) are only roughly indicated. The convex intruder where Sm.icro{l3,v) becomes multi-valued as function 
oi V > vp^ and /3 > /3p„ is indicated by shadowing. This corresponds to the first order transition. At the 
bottom the projection of the entropy surface onto the {/?, u} plane is shown as contour plot. The convex 
part (region of phase-separation) is hidden behind the dark "critical" line. 



Figure ( |10| ) explains what happens if one plots the entropy s vs. the "intensive" quantities 
f} — dS/dE and v — dS/dN as one would do for the grand-canonical ensemble: As there are 
several points Ei,Nk with identical Smicro{,l3,v) is a multivalued function of (3, v. Here the 
entropy surface Sraicrois-Tn) is folded onto itself. In the projection in fig.[To|, these points show up 
as a black critical line (dense region). Here this black line continues over the multi-critical point 
towards C indicating the direction to the critical point of the ordinary q = i Potts model at 
n = 1 (zero vacancies). Between and C the slopes 



ds 

d(3 



(31) 
(32) 
(33) 



are negative large but finite. 
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The information given by the projection would be all information which can be obtained from the 
conventional grand-canonical entropy s(T, /z, y), if we would have calculated it from the Laplace 
transform, eq.(||). The shaded region will be lost. 

The upper part of figure ( [To| ) shows Smicro{/3, in a three dimensional plot. The lines building 
the entropy surface are lines of equal f3. The images of the points A, D, B, C defined in fig.^ are 
roughly indicated. The back folded branches, the convex intruder of s(e,n) between the lines 
APmB and ADB, the region of phase separation, can here be seen from the side (shadowed). 
It is jumped over in eq. (||) and gets consequently lost in Z(r, /i). This demonstrates the far 
more detailed insight into phase transitions and critical phenomena obtainable by micro-canonical 
thermo-statistics which is not accessible by the canonical treatment, c.f. the similar arguments of 
Gibbs 

In the next two figures the cross-section through s(/3, ly) at constant (3 is shown in figure ( |Tl| ) 
below the multi-critical point f3m — 1-48 and in figure (^2|) above it. The latter clearly shows the 
back-bending of s(/3, i'). 




3 3.2 
FIG. 11. Plot of the entropy s(/3 = 1.325, v) 



1 . 1 




FIG. 12. Plot of the entropy = 1.7, v) 



VII. CONVEX ENTROPY — VIOLATION OF THE SECOND LAW ? 

At this point it is worth-wile to spend some words on a popular misunderstanding connected 
with the eventual convexity of the entropy as function of "extensive" quantities like the energy: 
The convex parts of S{E, N) violate van Hove's concavity condition [p^pl|]. 



18 



One may believe that this is also in contradiction to the second law of thermodynamics: At 
a convex region of S{E,N) a split of the system into two pieces with entropies Si{Ei, Ni) and 
S2iE2,N2) would have Si{Ei,Ni) + S'2(£^2, ^^2) > S{Ei + E2, Ni + N2). So the system seems to 
gain entropy by splitting. 

This, however, is an error. The Boltzmann entropy as defined in the section II is already the 
logarithm of the sum over all possible configurations of the system at the given energy. The split 
ones are a subset of these. Their partial phase space Wsput is of course < the total W. The entropy 
Sspiit — ln{Wsput) is < the total entropy. Evidently, the split system looses some surface entropy 
Ssurf at the separation boundary due to additional correlations imposed on the particles at the 
boundary, see the discussion in section fu[ The entropy after split is consequently: 

Sspiit = Si{Ei,Ni) + S2{E2,N2) — Ssurf 

<S{Ei+E2,Ni+N2), (34) 

It is a typical finite size effect. Ssurf /V vanishes in the limit V 00 for interactions with finite 
range. The entropy is non-extensive for finite systems but becomes extensive in the limit, and van 
Hove's theorem jl^ is fulfilled. This is of course only under the condition that limy^oo Ssurf /V = 
0. 

In general this is of course a trivial conclusion: An additional constraint like an artificial cut of 
the system can only reduce phase space and entropy. The Second Law is automatically satisfied 
in the Boltzmann formalism whether S is concave or not, whether S is "extensive" or not. 

A positive (wrong) curvature introduces problems to the geometrical interpretation of thermody- 
namics as formulated by Weinhold |^,^ which relies on the non-convexity of S{E, N). Weinhold 
introduces a metric like 

d^S 

where we identify : 
=E 
= N. 

The thermodynamic distance is defined as : 



Aa^b = ^J[X^{a) - X^{b)]g,k[X''{a) - Xfe(6)]. (36) 

Evidently, a negative metric gik is here not allowed. Of course Weinhold's theory does not apply 
to finite systems with phase transitions. 



VIII. CONCLUSION 



Micro-canonical thermo-statistics describes how the entropy s(e,n) as defined entirely in me- 
chanical terms by Boltzmann depends on the conserved "extensive" variables: energy e, particle 
number n, angular momentum L etc. It is well defined for finite systems without invoking the 
thermodynamic limit. Thus in contrast to the conventional theory, we can study phase transitions 
also in "small" systems or other non-extensive systems. In this simulation we could classify phase 
transitions in a "small" system by the topological properties of the determinant of curvatures 
(i(e,7i), eq.([l^) of the micro-canonical entropy-surface s(e,n): 

• A single stable phase by c?(e, n) > 0. 

• A transition of first order with phase separation by d{e, n) < 0. The depth of the intruder is 
a measure of the intra-phase surface tension |34|j2^ ]. This region is bounded by a line with 
d(e, n) = 0. On this line Pm is a critical end-point where additionally vi • Vd = in the 
direction of the eigenvector of d{e,n) with the largest eigenvalue Ai. 

• There, the transition is continuous ("second order") with vanishing surface tension, and no 
convex intruder in s{e,n). Here two neighboring phases become indistinguishable, because 
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there are no interfaces. However, we found a further line {PmC, critical) with Vi • Vd = 
which does not border a region of negative d{e, n). Presumably d(e, n) should be also. This 
needs further tests in other systems. It may also be that these lines signalize transitions of 
first order in other, but hidden conserved degrees of freedom. 

• Finally a multi-critical point P„i where more than two phases become indistinguishable by 
the branching of several lines with d = or with vi • Vd = to give a flat region with 
additionally = 0. 

Our classification of phase transitions by the topological structure of the micro-canonical Boltz- 
mann entropy s{e,n) is close to the natural experimental way to identify phase transitions of 
first order by the inhomogeneities of phase separation boundaries. This is possible because the 
micro-canonical ensemble does not suppress inhomogeneities in contrast to the grand-canonical 
one, as was emphasized already by Gibbs |2^. Intra-phasc boundaries are reflected in "small" 
systems by the convex intruder in the entropy surface. With this extension of the definition of 
phase transitions to "small" systems there are remarkable similarities with the transitions of the 
bulk. Moreover, this definition agrees with the conventional definition in the thermodynamic limit 
(of course, in the thermodynamic limit the largest curvature Ai approaches from above at phase 
transitions of first order) . The region of phase separation remains inaccessible in the conventional 
grand-canonical ensemble. 

We believe, however, that the various kind of transitions discussed here have their immediate 
meaning in "small" and non-extensive systems independently whether they are the same in the 
thermodynamic limit (if this then exist) or not. For systems like the Potts model that have a 
thermodynamic limit it might well be possible that the character of the transition changes towards 
larger system size. 

The great conceptual clarity of micro-canonical thermo-statistics compared to the grand- 
canonical one is clearly demonstrated. Not only that, we showed that the micro-canonical statistics 
gives more information about the thermodynamic behaviour and more insight into the mechanism 
of phase transitions than the canonical ensemble: About half of the whole {E,N} space, the in- 
truder of S{E,N) or the region between the ground state and the line APmB in figure (jj), gets 
lost in conventional grand-canonical thermodynamics. Without any doubts this contains the most 
sophisticated physics of this system. We emphasized this point already in |^ there, however, 
with still limited precision. Due to our refined simulation method this could be demonstrated 
here with uniformly good precision in the whole {E^N} plane. Finally, we should mention that 
micro-canonical thermo-statistics allowed us to compute phase transitions and especially the sur- 
face tension in realistic systems like small metal clusters [0. Our finding clearly disproves the 
pessimistic judgement by Schrodinger |]35| who thought that Boltzmann's entropy is only usefuU for 
gases. A recent application of micro-canonical thermo-statistics to thermodynamically unstable, 
collapsing systems under high angular momentum is found at 

Acknowledgment: D.H.E.G. thanks M.E.Fisher for the suggestion to study the Potts-3 model 
and to test how the multicritical point is described micro-canonically. We thank H. Jaqaman for 
critical reading. We are gratefully to the DFG for financial support. 



[1] T.L. Hill. Thermodynamics of Small Systems I+II. W.A.Benjamin Inc., New York, Amsterdam, 1964. 

[2] R.K. Pathria. Satisitical Mechanics. Pergamon Press, Oxford, New York, Toronto, 1972. 

[3] Elliott H. Lieb and J. Yngvaso n. The physics and mathematics of the second law of thermodynamics. 

Physics Report, 310:1-96, 1999, |cond-mat/970820c| . 
[4] L. Boltzmann. Uber die Beziehung eines algemeinen mechanischen Satzes zum Hauptsatz der 

Warmelehre. Sitzungsbericht der Akadamie der Wissenschaften, Wien, 11:67-73, 1877. 
[5] A. Einstein. Zur allgemeinen molekularen Theorie der Warme. Annalen der Physik, 14:354-362, 1904. 
[6] C. Tsallis. Possible generalization of Boltzmann-Gibbs statistics. J.Stat.Phys, 52:479, 1988. 



20 



G.A. Raggio. Equ ivalence of two thermostatistical formalisms based on the Havrda&Charvat-Daroczy- 



Tsallis entropies, http://xxx.lanl.gov/abs/cond-mat/9908207, pages 1-4, 1999 



J.W. Gibbs. Elementary Principles in Statistical Physics, volume II of The Collected Works of 
J. Willard Gibbs. Yale University Press 1902, also Longmans, Green and Co, NY, 1928. 
A. Einstein. Eine Theorie der Grundlagen der Thermodynamik. Annalen der Physik, 11:170-187, 1903. 
P. Ehrenfest and T. Ehrenfest. Begriffliche Grundlagen der statistischen Auffassung in der Mechanik, 
volume IV. Enzycl. d. Mathem. Wissenschaften, 1912. 

P. Ehrenfest and T. Ehrenfest. The Gonceptual Foundation of the Statistical Approach in Mechanics. 

Cornell University Press, Ithaca NY, 1959. 

Chapter VII of fsl, footnote on page 75. 

T. D. Lee and C. N. Yang. Phys. Rev., 87:410, 1952. 

D.H.E. Gross. Statistical decay of very hot nuclei, the production of large clusters. Rep. Progr. Phys., 
53:605-658, 1990. 

M. Schmidt, R. Kusche, W. Kornmiiller, B. von Issendorff, and H. Haberland. Experimental determi- 
nation of the melting point and specific heat for a free cluster of 139 sodium atoms, preprint Freiburg, 
1997. 

H. Haberland. Backbending and first order transition in small sodium clusters, private communication, 
1999. 

D.H.E. Gross and M.E. Madjet. Fragmentation phase transition in atomic clusters IV — the relation 
of the fragmentation phase transition to the bulk liquid-gas transition. Z. Physik B, 104:541-551 and 



http://xxx.lanl.gov/abs/cond-mat/9707100l , 1997. 

L. Boltzmann. Uber die Begriindung einer kinetischen Gastheorie auf anziehende Krafte allein. Wiener 
Berichte, 89:714, 1884. 

L. van Hove. Quelques proprietes generales de I'integrale de configuration d' un systeme de particules 
avec interaction. Physica, 15:951, 1949. 

J. P. Straley and M.E. Fisher. Three-state potts model and anomalous tricritical points. J. Phys. A. 
Math., Nucl. Gen, 6:1310-1326, 1973. 

J.W. Gibbs. Graphical methods in the thermodynamics of fluids, volume I of The scientific papers of 
J. Willard Gibbs. Longmans, Green, NY, 1906. 

D.H.E. Gross, A. Ecker, and X.Z. Zhang. Microcanonical thermo dynamics of first order phase 



transitions st udied in the potts model. Ann. Physik, 5:446-452, and http://xxx.lanl.gov/abs/cond- 
mat/9607150| 1996. 



F.Y. Wu. The potts model. Rev.Mod.Phys, 54:235, 1982. 

R.B. Griffith. Rigorous Results and Theorems, volume 1 of Phase Transitions and Critical Phenomena, 

chapter 1, pages 7-105. Academic Press, London, New York, (1972). 

D.A. Lavis and G.M. Bell. Statisitcal Mechanics of Lattice Systems. Springer, 1999. 

R.J. Baxter. J. Phys., C6:L445, 1973. 

A. Hiiller. Finite size scaling at first order phase transitions ? Z.Phys.B, 95:63-66, 1994. 

D.H.E. Gross and E. Votyakov. Phase transitions in finite svstems = topo logical peculiarities of the 



microcanonical entropy surface, http://xxx.lanl.gov/abs/cond-mat/9904 075 , page 4, 1999) 



E. Votyakov, H. Jaqaman, and D.H.E. Gross. Total microcanonical entropy surface s(e,n), continuous 
and first order phase transitions, critical lines and multicritical points for various finite q-state potts 
lattice gases. HMI to be published, 1999). 

I.D. Lawrie and S.Sarbach. Theory of tricritical points. In C. Domb and J.L. Lebowitz, editors. Phase 
Transitions and Critical Phenomena, number 9, chapter 1, pages 2-155. Academic Press, New York, 
1984. 

T.L. HiU. J. Phys. Chem., 23:812, 1955. 

F. Weinhold. Metric geometry of equilibrium thermodynamics. J. Chem. Phys., 
63:2479,2484,2488,2496, 1975. 

Frank Weinhold. Geometrical aspects of equilibrium thermodynamics. In H. Eyring and D. Hender- 
son, editors. Theoretical Chemistry, Advances and Perspectives, number 3, chapter 2, pages 16-52. 
Academic Press, New York, 1978. 

K. Binder. Monte Carlo calculation of the surface tension for two- and three-dimensional lattice-gas 
models. Phys. Rev. A, 25:1699-1709, 1982. 

E. Schrodinger. Statistical Thermodynamics, a Course of Seminar Lectures, delivered m January- 
March 1944 the School of Theoretical Physics. Cambridge University Press, London, 1946. 
V. Laliena. The efi'ec t of angular momentum conservation in the ph ase transitions of collapsing systems. 



PRE, 59:4786-4794, http://xxx.lanl.gov/abs/cond-mat/9806241 



1999. 



21 



